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Coherence, multiple scattering and the interplay between soft and hard processes are 
discussed. These physics phenomena are essential for understanding the nuclear depen- 
dences of rapidity density and pt spectra in high-energy heavy-ion collisions. The RHIC 
data have shown the onset of hard processes and indications of high spectra suppression 
due to parton energy loss. Within the pQCD parton model, the combination of azimuthal 
anisotropy (i^) and hadron spectra suppression at large px can help one to determine the 
initial gluon density in heavy-ion collisions at RHIC. 

1. Introduction 

Monte Carlo models have played an important and irreplaceable role in both low and 
high-energy physics. First of all, they are needed to simulate the physics and detector 
capabilities for any proposed experiment. They are necessary to study the efficiency, 
acceptance corrections and background in the analyses of experimental data. They are 
also very useful tools for theorists to test and estimate the consequences of new physics 
ideas and phenomena. This is especially true for high-energy heavy-ion collisions. 

In order to study the properties of the quark-gluon plasma (QGP) produced in the 
early stage of heavy-ion collisions and to search for evidence of the QCD phase transition, 
one needs to understand the whole evolution history of the collisions. Since the reaction 
dynamics is very complex and there is not a simple standard analytic model, one has to 
rely on Monte Carlo models to incorporate many aspects of strong interactions in the 
simulation of heavy-ion collisions. The models and parameters therein can be constrained 
by well-tested theories and experimental data especially in pp and pA collisions where 
QGP is not expected to form. They provide a baseline calculation of physical observables 
in the absence of new physics due to the formation of QGP. However, large uncertainties 
exist in extrapolating to AA collisions. Comparisons of results produced by different 
Monte Carlo models or by varying the model parameters are useful to provide a measure 
of the extrapolation uncertainties. The latest RHIC data |]-§J have proven to provide 
good constraints on Monte Carlo models that in turn help us to understand the reaction 
dynamics and the initial conditions. Finally, Monte Carlo models also serve as theoretical 
laboratories to test proposed signals and probes of QGP such as jet quenching 0. 
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There are many Monte Carlo models for heavy-ion collisions, reflecting the complexity 
of the problems. But in general they can be divided into three categories. Hadron and 
string based models, e.g., FRITIOF §, DPM |], VENUS [|, RQMD @, ARC 0, ART 
p"T|j , URQMD |nj, LUCIAE [|13|] and others, consider hadrons and Lund strings ]T3[] as 
the effective degrees of freedom in strong interactions. Depending on the center-of-mass 
energy of a two-body collision, particles are produced through resonances and strings 
which hadronize according to the Lund string fragmentation model |14|]. These models 



are adequate for particle production and rescattering in hadronic and nuclear collisions 
at energies around or below the CERN-SPS, in which partonic degrees of freedom are not 
yet important. At high energies (y'i > 50 GeV), hard or semi-hard partonic scatterings 
become important and contribute significantly to the particle production and dominate 
the high px hadron spectra even at the CERN-SPS energy. In this case, pQCD-inspired 
models, e.g., HIJING |TBJ, VNI [|T5| NEXUS g7| and AMPT [T8fl, are more relevant. The 



third class of models, such as LEXUS fL9|l , pQCD parton model pO|,^TJ, parton cascade 
models P^fl and hydrodynamic models 



], are not strictly event generators. But they are 
useful to help us to understand particle production and evolution in heavy-ion collisions, 
many aspects of which are difficult to simulate in an event generator. 

I will not give a review of these models [24| here. Instead, I will concentrate on the 



parton based models and discuss the essential physics behind them. In particular, I will 
discuss the interplay between soft and hard processes and their different behaviors in 
multiple parton scattering processes. I will demonstrate the difficulties in incorporating 
quantum interferences in Monte Carlo models and how one can use mission specific models 
such as the pQCD parton model to overcome these difficulties. 

2. Soft versus hard processes 



|] and heavy-ion 

was realized, two parton-based Monte Carlo models pT5|,p!6| were developed 



Not long after the importance of minijets in both high-energy pp{p) 
|26 



collisions 

to incorporate hard processes in high-energy heavy-ion collisions. The jet production cross 
section in these hard processes can be described very well by the pQCD parton model 
when the transverse momentum transfer involved is very large. However, the pQCD 
parton model calculation of jet cross sections diverges when — > and the processes 
become non-perturbative. No theoretical model exists that can describe both the large 
and small pt processes. In order to combine hard processes and non-perturbative soft 
processes in a Monte Carlo model, a cut-off scale po in transverse momentum transfer of 
parton scattering is introduced that phenomenologically separates soft and hard processes. 
Assuming eikonalization of these processes, one can calculate the pp(p) cross section in 
this two-component model 
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where T NN (b) is the nucleon-nucleon overlap function, cr so f t is a parameter to represent 
inclusive soft parton cross section and <Tj e t is calculated from pQCD parton model, 
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The two parameters, po and <Tj et in this two- component model can be fixed by the experi- 
mental data on the cross sections of pp(p) collisions and the corresponding dN^/drj. This 
is the basic principle behind HIJING @, VNI |L6|] and the recent NEXUS 0. In this 



two-component model, the rapidity density of charged multiplicity can be expressed as 

In HIJING, soft particle production is modeled by a string model whose contribution 
to dN^/drj is almost constant. The contribution from minijet production, on the other 
hand, increases with energy due to the increase of the inclusive jet cross section. 

To extend to nuclear collisions, one has to consider the nuclear dependences of these 
two components. For the soft part, the interaction is coherent and its contribution to the 
multiplicity should be proportional to the number of wounded nucleons. On the other 
hand, minijet production is approximately incoherent and proportional to the number of 
binary collisions, modulo the nuclear shadowing of parton distributions and jet quenching. 
The total charged hadronic rapidity density will acquire the form 

1 dNtf _ ,v ,,v (AW) m 
<AWt>/2 d V Wsoft + Whard (iV part >/2 ' [ ) 

where a^(s) is the averaged inclusive jet cross section per NN collision in AA, which 
is generally smaller than a™ due to nuclear shadowing of parton distributions. The 
energy dependence of the minijet cross section a^(s) in the hadronic multiplicity is then 
amplified by (Aginary) / (Apart) ~ A 1 / 2, in heavy-ion collisions relative to pp(p) by the binary 
nature of semi-hard processes. 

Shown in Fig. [I] are HIJING results of dN c h/ drj per pair of participant nucleons in pp(pp) 
and central Au + Au collisions at different colliding energies as compared to experimental 
data. The rise of the multiplicity per participant relative to SPS is consistent with the 
predicted binary scaling of the hard component from pp to AA via Eq. @ without jet 
quenching. Final state interaction such as jet quenching will enhance the A dependence 
of dN ch /dr] especially at high energies. 

Also shown in Fig. [TJ is the result of a final-state saturation model by EKRT |]29|| . In 
this model, the pQCD growth of low pr gluons is cut off below a saturation scale which 
is determined by the local saturation requirement, dN g /dyd 2 s = p 2 sat /^- Here dN g /dyd 2 s 
is the gluon rapidity density per unit transverse area. The resultant saturation scale 
increases with energy and (N part ). Assuming direct proportionality between parton and 
the final hadron number, the energy and nuclear dependence of the total hadronic rapidity 
density in central A + A collisions can then be estimated. As shown, the EKRT result is 
remarkably close to the HIJING results and is also consistent with the RHIC data. 

However, a critical difference between these models is the centrality dependence of 
dN C h/drj(N pai t). In the HIJING two-component model of minijet production the multi- 
plicity per participant increases with A according to Eq.(|j) while it actually decreases in 
the EKRT model with increasing A due to the saturation requirement. In order to distin- 
guish the fixed-scale and saturation models of entropy production, the study of centrality 
dependence of the hadronic multiplicity per participant was proposed . 
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Figure 1. The RHIC data 0-0 for central 
Au+Au are compared to pp and pp data 
and the NA49 Pb + Pb data 
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28 . HI- 



JING1.35 (solid) with (upper) and with- 
out jet quenching (lower) and EKRT (dot- 
dashed) predictions |29| are also shown. 



Figure 2. The RHIC data @,f§ on central- 
ity dependence of dN ch /drj/0.5(N pSbTt ) are 
compared with HIJING and EKRT |29 
predictions. The two-component model in 
Eq. (|]) with (solid) and without (dashed) 
b- dependence of shadowing are also shown. 



Shown in Fig. § are the dN c h/dr] per participant pair as functions of (N paxt ). The 
EKRT model [2S] has a weakly decreasing or constant dependence on centrality in semi- 
peripheral to central collisions. The HIJING result, on the other hand, increases mono- 
tonically with the number of binary collisions per participant (iV binai . y ) / (N pait ) in an in- 
tuitive way as given in Eq.([|). The RHIC data clearly favor the the latter case. Naively, 
(iVbinary)/(-^part) ~ (^ pa rt) 1//3 - The deviation from such a simple dependence in HIJING 
calculation might be due to a combined effect of jet quenching and dilute edges in nuclear 
distributions used in HIJING. Also shown as dashed line is the simple two-component 
model in Eq. (Q) which agrees with the HIJING result. The parameters, (n) so f t « 1.6 and 
(^)hard ~ 2.2 are determined from a fit to the measured energy dependence of dn^/dr] in 
pp(p) collisions. Another phenomenological two-component model by Kharzeev and Nardi 
[ |32]| , shown as the dotted line, assumed dn^/drj = 2.25 when fixing the parameters. It 
coincides exactly with the dashed line if dn^/dr] = 2.0 is used. 

The jet cross section crdf used in the two-component model (dashed line) is calcu- 
lated with nuclear shadowing of parton distributions averaged over mini-biased events. 
Shadowing reduces the effective jet cross section and thus minijet contribution to charged 
multiplicity. Without shadowing, the two-component model will over-predict the charged 
multiplicity in central Au + Au collisions. However, parton shadowing should depend 
on the impact-parameter of nuclear collisions. If an impact-parameter dependent parton 
shadowing as described in HIJING |TJj is used, the effective jet cross section a££ should 
also depend on centrality. Using this 6-dependent cross section in Eq. (U), one gets the 
solid line in Fig. 0. The RHIC data certainly favor this scenario. Such shadowing effect 
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is similar to the initial state saturation [30]. A phenomenological implementation of the 



initial state saturation is shown to give a similar result as the two-component model. 
Another independent constraint is the rapidity dependence of dN ch /dr]/ (N part ). PHOBOS 
preliminary [H] data are consistent with HIJING model. 



3. Interplay between soft and hard processes 

The centrality dependences of soft and hard processes not only manifest themselves in 
the total charged multiplicity but also in the hadron spectra. Incoherent hard parton 
scattering should dominate the hadron spectra at high while coherent soft interactions 
contribute mainly to the low pt region. This leads to some nontrivial nuclear dependence 
of hadron p T spectra. A good Monte Carlo model should be able to simulate this nuclear 
dependence. 

In a schematic Glauber model of multiple parton scattering, one can calculate the 



produced parton spectra 33] in p + A collisions: 



Edo 3 



NA 



A d 3 p 



E 



da 3 

aa NN 

d 3 p 



QA I' v r dxi . . 

^T^L / —fi/N{Xi) 
J A j 



lQlTTn 



E 



d 3 p k 
Eh 



u iN u kN 



(a iN + a jN )hl 



iN 



(5) 



where h\ N is the differential cross section for parton-nucleon scattering i + N — > j+X, and 
/i/jv(^t) is t ne parton distribution in a nucleon. The effective parton-nucleon total cross 
section is defined as <JiN(Pi) = \ Hj I -f^^iNiPiiPj) an d the differential nucleon-nucleon 
cross section for parton production as Eda J NN /d 3 p = Y,i f fi/N(xi)h J iN (pi, p) . We have 
taken a nucleus to be a hard sphere of radius ta = r^A 1 / 3 . 

In the above equation, the first term corresponds to the incoherent sum of single par- 
ton scatterings in a nucleus. The second term gives the nuclear modification of the 
parton spectrum due to multiple parton scattering inside a nucleus. This term con- 
tains contributions from both the double parton scattering and the negative absorptive 
correction. For a schematic study, let us assume that all partons are identical and 
the differential parton-nucleon cross sections have a simple regularized power-law form 
h\ N = h{p T ) = C/(pt + Po) n . The parameter p can be considered as the scale separating 
soft and hard processes. One can then calculate the nuclear modification factor of the 
parton spectra, 



Ra = 



EdaNA/d 3 p 
AEdaNN I d 3 p' 



(6) 



which is shown in Fig. ^ for three different values of as a function of Pt/Po- Relative 
to the additive model of incoherent hard scattering, the spectra are enhanced at large pt 
(Ra > 1) for hard processes due to multiple scattering. For soft processes at low px, the 
absorptive correction is large and the second term in Eq. (JD becomes negative, leading 
to the suppression of hadron spectra (Ra < 1)- The transition between suppression and 
enhancement occurs at around the scale po that separates soft and hard processes. If we 
analyze the hadron spectra in pp(p) collisions as shown in Fig. [|, one can indeed see the 
two underlying components: one in power-law form that dominates spectra at high px 
and another in exponential form for low pt hadrons. The transition between these two 
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Figure 3. The nuclear modification factor in Figure 4. Experimental data |27j on hadron 

a schematic model of multiple parton scat- spectra in pp(p) collisions are compared to 

tering with a simple form of parton-nucleon two-component (power-law and exponen- 

cross section l/(j4 + Po) n . tial) parameterization. 



components occurs at around po 
collisions 



20,33 



1.5 — 2 GeV/c. Analyses of the hadron spectra in pA 
in terms of the nuclear modification factor indeed show such a feature 
in Ra{pt) and the transition scale is about p « 1.5 GeV/c. These same features from 
initial multiple parton scattering should also appear in heavy-ion collisions if there is no 
final state interaction. Shown in Fig. |5] are the experimental data of central Pb + Pb 
collisions at the CERN-SPS. These data have almost exactly the same behavior as in 
pA collisions at the same energy. The soft-hard transition scale is also around po ~ 1.5 



GeV/c. However, there seems to be no evidence of jet quenching at SPS |20] which would 
suppress the high p T spectra. 

As demonstrated in the schematic model, the unique features of the nuclear depen- 
dence of the hadron p? spectra are consequences of the intricate quantum interference 
between multiple parton scattering and absorptive processes. At low p?, the absorptive 
processes are dominant and they reduce the contribution of independent parton scatter- 
ing. Consequently, the hadron production at low pt appears to be coherent. The inclusive 
differential cross sections are only proportional to the surface area of the colliding nuclei or 
the number of participant nucleons. At large pt, the destructive interference between mul- 
tiple parton scattering and absorptive processes is almost complete. Hadron production 
processes at large p? become incoherent. The corresponding inclusive cross sections are 
proportional to the volume of the colliding nuclei or the number of binary nucleon-nucleon 
collisions. It is difficult to simulate this quantum interference in a Monte Carlo model, 
especially in the transition region where both processes are important. In HIJING, the 
coherent soft processes are modeled by Lund-like strings whose number is proportional to 
the number of participant nucleons. The hard processes are treated as incoherent parton 
scatterings (modulo nuclear shadowing of parton distributions) that are proportional to 
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Au+Au(b=0) Vs=130GeV 




Figure 5. The nuclear modification factor Figure 6. The nuclear modification factor 

for hadron spectra in central Pb + Pb col- for hadron spectra in central Au + Au colli- 

lisions at the CERN-SPS. The solid line is sions at RHIC with different values of par- 

pQCD parton model calculation ton energy loss. 



the number of binary collisions. Because of parton fragmentation, hard processes still 
contribute to hadron production at low px- However, since a px cut-off at p ~ 2 GeV/c 
is imposed in HIJING, hard contributions to hadron production become smaller at lower 
Pt- Soft processes also contribute to hadron production with an exponential-like spec- 
trum. Therefore, the resultant hadron spectra in HIJING have a smooth transition from 
soft to hard processes as pt is increased f| . However, the string fragmentation used for 
jets tend to suppress the pt broadening due to multiple parton scattering. So HIJING 
cannot reproduce the Cronin enhancement of hadron spectra at intermediate px- 



4. Parton model and high p T spectra 

The transition from soft to hard processes is only one example of quantum interferences 
one has to consider in a Monte Carlo model for high-energy heavy-ion collisions. Other 
effects include LPM interference effect for a propagation of an energetic parton inside a 
medium, coherent multiple scattering in a dense medium, and effects of chiral symmetry 
restoration. Implementing these phenomena correctly in a Monte Carlo model that can 
simulate every stage of heavy-ion collisions is very difficult, if not impossible. We can, 
however, resort to some mission-specific models to help us to understand the physics 
of some specific aspect of heavy-ion collisions. For example, hydrodynamic models |23[ 
can help us to understand the collective phenomena in situations where cascade models 
become invalid because of coherent scattering in a dense medium. Parton cascade models 



22] can assist us in understanding the process of parton thermalization. In the following, 
I will discuss the pQCD parton model P0| , pT| , which is very useful for the study of high px 
spectra from jets. They can be used as a probe of the dense partonic matter in heavy-ion 
collisions. 

For high px hadron spectra, one can essentially neglect the soft processes and only 
consider hadron production from hadronization of jets whose cross section can be de- 
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Figure 7. The preliminary 7r spectra in cen- 
tral Au + Au at RHIC || as compared to 
pp data |27| and parton model calculations. 
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Figure 8. The hydro+parton model calcu- 
lation of V2{pt) as compared to preliminary 
STAR data M. 



scribed by Eq. (0). Hadron spectra in pp, pA collisions can be described well j2(| within 
this parton model. Constrained by the existing pA data, nuclear modification of parton 
distributions and px broadening give about 10-30% increase in the pt spectra in central 
Au + Au collisions relative to pp without parton energy loss, as shown in Fig. 

Many recent theoretical studies |34| predict a large amount of energy loss by a prop- 
agating energetic parton in a dense partonic medium. In AA collisions, we model the 
effect of parton energy loss by the modification of parton fragmentation functions . 
This will lead to suppression of large pt hadrons for a non- vanishing parton energy loss 
relative to the null scenario. In terms of the nuclear modification factor, Rab{pt) will 
become smaller than 1 at large pt due to parton energy loss as shown in Fig. |6|. It is also 
shown that the pr dependence of Rab(pt) is sensitive to the energy dependence of the 
parton energy loss. For more accurate determination of the energy dependence of dE/dx 
one should measure the high px hadron suppression in direct 7-tagged events |35|| . 

Preliminary RHIC data have indeed shown evidence of suppression of large px 
hadron spectra for the first time. This is extraordinary, given that no suppression is seen 
at SPS f20fl . Shown in Fig. |7] are 7r spectra in central Au + Au collisions as compared 
to pQCD parton model calculation. The model result for pp collisions at yfs = 200 GeV 
agrees well with the UAl data. Both the calculation and data for Au + Au collisions 
are scaled by (N paxt ) = J d 2 bTAB(b) for the given centrality in order to compare to pp{p) 
results. The preliminary data clearly show suppression of high px hadron spectra relative 
to pp. The data are consistent with the parton model calculation with an average dE/dx ~ 
0.25 GeV/fm. I want to emphasize that this is only a phenomenological value averaged 
over the entire evolution history of the dense matter. It does not reveal the distance and 
density dependence of the dE/dx. Since the initial parton density decreases very fast 



9 



due to longitudinal expansion, such a small average dE/dx still implies very large parton 
energy loss in the very early stage of the dense system. Future analysis of the energy and 
centrality dependence |36| of the high p T hadron suppression is important to understand 
the consistency between CERN-SPS and RHIC results and search for the critical initial 
parton density at which jet quenching becomes significant in heavy-ion collisions. 

In non-central AA collisions, the total parton propagation length should depend on the 
azimuthal direction. Therefore, parton energy loss can also cause azimuthal anisotropy of 



hadron spectra at large px in non-central AA collisions |j36| . Such azimuthal anisotropy 
can be calculated in the same parton model with a given parton energy loss. Shown in 
Fig. H is the calculated |36j v 2 {j>t) as a function of pr as compared to STAR prelimi- 
nary data Q. The hadron spectra at low px < 2 GeV/c are assumed to be given by 
the hydrodynamic model and v 2 increases linearly with p?. However, at pt > 2 GeV/c, 
hydrodynamics will fail and v 2 should be determined by the dynamics of parton propa- 
gation in the dense medium. The parton model gives very different pt dependence. The 
point where v 2 (pt) deviates significantly from the hydro model signals the onset of the 
contribution of hard processes. As shown in the figure, v 2 at large pt is very sensitive to 
dE/dx or the initial parton density. It is an alternative measurement of parton energy 
loss or the initial parton density. Since suppression of hadron spectra and non- vanishing 
v 2 at large pt are two consequences of parton energy loss, one should be able to explain 
both experimental measurements and their centrality dependence with the same given 
parton initial density. 



5. Summary 

I have discussed the importance of multiple scattering, coherence and the interplay 
between soft and hard processes in the modeling of high-energy heavy-ion collisions. Both 
the nuclear dependences of rapidity density and pt spectra have shown the onset of hard 
processes. The preliminary RHIC data show evidence of high p? hadron suppression due 
to parton energy loss. I argued that the combination of v 2 and hadron spectra suppression 
will help us to determine the initial gluon density in heavy-ion collisions at RHIC. 
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